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ABSTRACT 

Cosmic ray (CR) acceleration at the shock created by the expanding cocoons 
around active galactic nuclei (AGNs) is studied. It is shown that above the 
energy 10 18 eV the overall energy spectrum of CRs, produced during the AGN 
evolution and released in the intergalactic space, has the form N oc e -7 , with 
7 « 2.6, which extends up to e max ~ 10 20 eV. It is concluded that cocoons shocks 
have to be considered as a main source of extragalactic CRs, which together with 
Galactic supernova remnants provide the observed CR spectrum. 

Subject headings: cosmic rays - - acceleration of particles — shock waves - 
supernova remnants — galaxies: active 



1. Introduction 

The overall origin of cosmic rays (CRs) is still an unresolved problem in astrophysics. 
The understanding of CR origin requires the determination of the astrophysical objects, that 
are the CR sources, and of the appropriate acceleration processes, that form the CR spectrum 
in these objects. In this regard, during the last several years considerable progress has been 
achieved in this field, both experimentally and theoretically . Recently the steepe ning of the 



CR spectrum aboye 3 x 1 eV was established in the HiRes (IBergman et al.ll2007l ) and Auger 



(lYamamoto et al.l 120071 ) experiments. It presumably corresponds to the so-called Greizen- 
Zatsepin-Kuzmin (GZK) cutoff, caused by CR energy losses in their interactions with the 
microwave background radiation. This is evidence that the highest ener gy part of the CR 



spect rum is of extragalactic origin. It was also recently demonstrated (IBerezhko fc Yolk 



20071 ) that the CRs with energies up to e ~ 10 17 eV can be produced in supernova remnants 



(SNRs). 

The main reason why SNRs are usually considered as the CR source is a simple argument 
about the energy required to sustain the Galactic CR population against loss by escape, 
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nuclear interactions and ionization energy loss. Supernovae have enough power to drive the 
CR acceleration. The high velocity ejecta produced in the supernova explosion interacts 
with the ambient medium to produce a strong blast wave, which may accelerate a small 
suprathermal fraction of the ambient plasma to high energies. The only theory of particle 
acceleration which is sufficiently we ll developed and specific to allow quantitative model 



calculations is diffusive acceleration (IKrvmskv 



1977 



shock associated with SNRs (e.g. see lBerezhko fc Krymskylll988l ; iBerezhkol 12005k 120081 . for 



Bell 19 78T) applied to the strong o uter 



review) . Considerable efforts have been made during the last years to empirically confirm the 
theoretical expectation that the main part of CRs indeed originates in SNRs. Theoretically, 
progress in the solution of this problem has been due to the de velopment of a kinetic nonlinear 
theory of diffusive shock acceleration (IBerezhko et al.lll996l ). This theory includes all the 
most relevant physical factors, essential for the SNR evolution and CR acceleration in a 
SNR and it is able to make quantitative predictions of the expected properties of CRs 
pro duced in young SNRs and their nonthermal radiation. Applied to individual young SNRs 



; see iBerezhkol 120051 . 120081 . for reviews) this theory has successfully explained many observed 
SNR properties. Making use of the observed synchrotron emission spectrum from radio to 
X-ray frequencies it permits the determination of th e degree of magnetic field amplification, a 
process advocated earlier from plasma simulations (jLucek &: Belli l2000l : iBell fc Lucekll200ll ). 
which leads to the concentration of the highest-energy electrons in a very thin shell just 
behind the shock. Such observed filamentary structures have been observed i n har d X- 



rays (jVink fc Laming! 120031 ; lLong et al.l 120031 ; iBamba et al.l 120031 ; IParizot et al.l 120061 ) and 

are nowa days used as a second independent method to infer the magnitude of the amplified 
field 



sec 



Berezhko et al.ll2003l . for interpretation) 



It has been demonstrated, that magnetic field amplification leads to considerable in- 
crease of the maximal energy of CRs accelerated in SN Rs and that the observed CR energy 
spectrum can be well represented by two components (IBerezhko fc Volkl 120071 ) . The first 
one, J g (e), dominant up to 10 17 eV, consists of CRs produced in Galactic SNRs , whereas the 



secon d, J e9 (e), is produced in extragalactic sources. This so called dip scenario (lAloisio et al. 



20071 ) requires a relatively steep CR spectrum produced in extragalactic sources, J^ g oc e~ 7 
with 7 = 2.55 — 2.75 at energies e > 10 18 eV. It has been demonstrated that peculiari- 
ties in the CR spectr um - the so-called "dip" structure at e ~ 10 19 eV and a GZK cutoff 
at e > 3 x 10 19 eV (lAloisio et al.l 120071 ) - produced due to CR interaction with the cos- 
mic microwave background radiation on their way from the source to the observer are very 
consistent with the experiment. 

Below we consider CR acceleration at the shock created by the expa nding cocoon around 
active galactic nuclei (AGNs). Compared with previous considerations (INorman et al.lll995l ) 
magnetic field amplification is taken into account, which provides CR acceleration with an 
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appropriate spectrum up to the energy 10 20 eV. Therefore the cocoon shocks have to be 
considered as a main source of extragalactic CRs, which together with Galactic supernova 
remnants presumably provide the observed CR spectrum. 



Results and discussion 



Astrophysical objects, which are considered as potential extragalactic sources of ultra 
high energy CRs, should fulfill a number of requirements. First, they should have sufficiently 
high energy output, not less than L p ~ 2 x 10 — 3.5 x 10 erg M pc~ 3 yr -1 depending on 
the form of CR spectrum that they produce (IBerezinsky et al.l 120061 ) . 



Given the energy requirements, AGNs (e.g. IBerezinsky et al.l 120061 ) and gamma-ray 



bursts (GRBs) flMikrom & Usov! Il995l : IWaxmanl Il995 



Vietri 



19951 ) are considered as po- 



tentia l extragalactic sources of the ultra-high-energy CRs. However, as IBerezinsky et al. 



( 120061 ) argued, the energy output of GRBs has a serious problem if they are considered as 
the main source of the extragalactic CRs. 

The second requirement to the extragalactic CR sources is their ability to produce a 
power-law CR spectrum up to the maximal energy, which for protons should be at least 
as large as e max = 10 20 eV, which is well above the GZK cutoff. Particle acc eleration at 



relativistic shocks in AGNs (lAchterberg et al. 
( Milgrom fc Usov 1995 ; WaxmanT 



1995 



Vietri 



200d : iRachen & Biermannlll993l ) and GRBs 



19951 ) are discusse d as appropriate acc elera 



tion processes. However, simulations of CR shock acceleration (INiemiec et al.l 120061 ) have 
demonstrated a low efficiency in the case of relativistic shocks. 

The relativistic jet in AGNs is surrounded by the hot coc oon, which ex pands into the 
surrounding intergalactic medium with supersonic speed (e.g. iFerraril 119981 ) . Since the es- 



sential fraction of the jet energy is transformed into the internal energy of the background 
medium due to the outermost nonrelativistic shock, driven by the cocoon overpressure, it is 
natural to suppose that a good part of this energy is represented by effectively accelerated 
CRs, like what takes place in SNRs. Since CR acceleration by nonrelativistic shocks has 
been very well studied one can obtain reliable estimates of the expected spectrum of CRs. 
First of all the maximal energy of CRs, accelerated at the expande d shock of size R(t) and 
speed V = dR/dt, is determined by the expression (lBerezhkdll996l ) 



ZeBRV/c, 



where B is the upstream magnetic field, c is the speed of light, e is the proton charge and Z is 
the charge number of CR nuclei. This expression consistently takes into account the influence 
of all the most essential physical factors restricted the maximal CR energy, which are the 
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adiabatic cooling of CRs in the downstream region and the shock dece leration followed by 
the CR escape from the shock vicinity into outer space (IBerezhkol 119961 ) . 



Magne tic field near the shock front, as was established for all young SNRs (IBerezhko 



20051 . 120081 ) , is strongly amplified nonlinearly by CRs up to the level 



B 2 /(8ir) w 3 x lO^V^ 2 



(2) 



which is presumably also appropriate for extragalactic shocks. Here p is the external gas 



density. Note that the relation B oc V foll owing from Eq 
for the shock speeds V < 10 4 km s— 1 (IPelletier et a. 



amplified magnetic field, as was originally shown by lBell 



2) is indeed theoretically expected 
2006). Due to the nonresonantl y 



(120041 ) (see also IPelletier et al.ll2006h 



at higher shock speed V > 10 4 km s— 1 the expected relation is B oc V 3 ^ 2 , which leads to 
larger magnetic field strength compared with the value given by Eq.(2). Due to this fact CR 
maximal energy e max can be a few times larger compared with our estimate below. However 
this does not play a critical role because even our estimate gives e max ~ 10 20 eV, which is 
well above GZK cutoff. 



The outer cocoon shock expands with the speed (e.g. iFerraril 119981 ) 



(3) 



where Lj is the mechanical luminosity of the jet and VJi is the hot spot or the jet head speed. 
Substituting the magnetic field value following from Eq.(2), the shock speed (3) and the 
shock size R = 2Vt into the equation (1), we have 



10 20 Z 



u 



10 46 erg/s 



3/4 



N„ 



10~ 4 cm" 3 



-1/4 



10 10 cm/c 



-1/4 



1 kpc 



-1/2 



eV, 



(4) 



where Zj = Vht is the jet length. According to this expression, during the evolutionary epochs 
corresponding to the jet length from Zj = 1 kpc to Zj = 10 Mpc the cocoon shock produces 
a power-law proton spectrum which extends up to the maximal energy from e max = 10 18 eV 
to e max = 10 20 eV. Thus we conclude, that AGNs fulfill also the second requirement for 
extragalactic CR sources: protons, the main kind of ions in the space plasma, are effectively 
accelerated at least up to the energy e max = 10 20 eV. Heavier ions are expected to be 
accelerated to Z times higher maximal energy. However, since the GZK cutoff energy of 
nuclei e^zK < eozx for the atomic n umber A > 1 is lower then that for protons €qzk ~ 
5 x 10 19 eV (IBerezinsky et al.l 120061 ) , the predicted dependence e max oc Z is difficult to 
establish experimentally. 
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The most realistic possibility of finding the signature of the acceleration process of ultra- 
high-energy CRs is the experimental study of CR composition at energies e > 10 17 eV. If 
extragalactic CRs are indeed produced by large-scale nonrelativistic shocks their composition 
will have the following peculiarity. Inside and near the source, CR composition corresponds 
to the composition of the intergalactic gas with substantial enrichm ent by heavy elements , 



like with Galactic CR versus interstellar medium compositions (e.g. iBerezhko et al.lll996l ). 
If compositions of intergalactic and interstellar plasmas are not very different, then the 
expected source CR composition is characterised by the mean atomic number < A 5, like 
Galactic CR components. Compared with this source composition, the observed composition 
is significantly distorted at almost all energies due to propagation effects in the intergalactic 
medium and in the Galactic wind. Only at the energy e* ~ 2 x 10 18 eV are these compositions 
roughly the same because it is just below the lowest GZK cutoff energy e§| K ~ 4xl0 18 eV and 



the m odification of the proton spectrum at this energy is still not very large (IBerezinsky et al. 



20061 ). Since CR diffusion mobility for given energy e is lower for larger Z, heavier CRs 
undergo stronger modification (depression) during their propagation from the sources to 
the observer. Due to this effect composition of extragalactic CRs is expected to become 
progressively lighter with the decrease of their energy at e < e*. At higher energies, e > e*, 
due to the dependence of the GZK cutoff energy eczK oc A on the nuclear mass number 
A, CRs become progressively heavier so that the peak value o f < A > is expected a t the 
beginning of the GZK cutoff for iron, which is at e « 10 19 eV (IBerezinsky et al.ll2006l ). At 
such a high energies only protons and iron nuclei survive and the iron-to-proton ratio has 
a maximum value at the energy e ~ 10 19 eV, which provides a peak of CR mean mass 
very similar to the peak at e m 10 17 eV. To make a quantitative prediction of the expected 
ultra-high-energy CR composition, one has to perform a detailed consideration, which will 
be done in a subsequent paper. 

Note that the expression (4) is based on the amplified magnetic field value (2) and 
therefore does not depend o n the assumed value of intergalactic magnetic field in contrast to 



the earlier estimate of e max ( INorman et al.lll995l ). Note also that as is seen from the Eq.(4) 



maximal CR energy only slightly depends on external gas density and the jet head speed. 
The values of mechanical luminosity Lj and the jet length Zj, which determines the maximal 
CR energy, can be directly measured in the experiment. 

The third requirement to CR sources is related to the form of the CR spectrum. The 
form of the resulting CR spectrum produced during the whole evolution of the expanding 
shock is determined by three physical factors: (1) nonlinear shock modification due to the 
CR back-reaction, (2) adiabatic energy losses in the downstream region, and (3) diffusive CR 
escape from the shock vicinity into outer space. The existence of the CR escape phenomenon 
makes it possible to estimate the shape of the spectrum of the most energetic CRs produced 
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during the sour ce evolution. As was demon strated analytically (IBerezhko fc Krymskylll988 
Berezhkolll996l ) and confirmed numerically (IBerezhko et al.lll996l ). since the maximal energy 
e m ax(£) of CRs accelerated by the expanding shock during its most active phase decreases with 
time t due to shock deceleration, the most energetic particles with energy e ~ € ma _ x {t) undergo 
diffusive escape from the shock vicinity into the surrounding space. If such behavior happens 
during the time period from t\ to t 2 , then at any given epoch ti < t < t 2 all accelerated 
particles with energies e > e max (t) have already escaped and inside the source there are only 
particles with e < e max (t). 

Note that the particle escape phenomenon, has two important aspects. First, the en- 
ergy of escaped CRs determined for any evolutionary epoch by Eq.(l) does not depend 
any more on the adiabatic cooling inside the parent expanding source, like what happens 
with CRs of lower energies e < e max (t). Second, the spectrum of the escaped CRs can be 
essentially different compared with the canonical spectrum N cx e -2 , produced by strong 
nonrelati yistic shock. To estimate the spectrum of escaped CRs one can use the simple 
relation (IBerezhko fc Krymskyill988l ) 



N{e)ede oc pV 2 R 2 ~ /3 dR, 



(5) 



which determines at any given phase t the overall number of CRs N(e) with highest energy 
e ~ e max(£)- Due to the hard self-consistent spectrum of CRs, produced by the strong shock, 
the main co ntribution to the CR e nergy content is given by the particles with highest energy 
e ~ e max (t) (IBerezhko et al.lll996l ). Therefore their energy content scales as the shock ram 
pressure pV 2 times the shock volume, as it is given by Eq.(5). Factor R~@ describes the 
progressive slow decrease of CR acceleration efficiency due to the shock weakening. For the 
shock expansion law R oc t~ v relation (5) gives 



iV(e)oce- 7 with 7 = 1 + (2 - 0)/(2/i/ - 3). 



(6) 



In our case v — 1/2 this gives 7 = 3-/3. For the constant acceleration efficiency we have 
/3 = which gives 7 = 3. However, the shock deceleration is accompanied by the decrease 
of the acceleration efficiency. This effect is described by the amount of shock modification, 
which is characterized by the shock comp ression ratio, which f or the case of a strong shock 
depends on the shock speed as o oc V 3/8 (IBerezhko et alJll996r ). Since in our case V oc 1/R 
we have (3 = 3/8 that gives 7 ~ 2.6. Such a spect rum of extragalactic C Rs at energies 
e > 10 18 eV corresponds very well to the experiment (IBerezinsky et al.l 120061 ) . 



Particles with energies e < e max (£) at any given epoch are contained inside the source 
and have a spectrum close to iV oc e~ 2 . At a very late epoch when the outer shock becomes 
weak these particles will leave the source. Note however that this is already not important for 
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the Galaxy, because the contributio n of extragalactic CRs is expected to be low at energies 
e < 10 18 eV feerezinskv et al.ll2006h . 



Note that based on the data collected at the Auger experiment, a correlation between 
the arrival direct ions of CRs with ener gy above 6 x 10 19 eV and the position of nearby AGNs 
has been found ([Abraham et al.l 120071 ). which strongly supports AGNs as prime candidates 
for the source of ultra-high-energy CRs. 



3. Summary 

Magnetic field amplification due to shock-accelerated CRs provides acceleration of CRs 
up to the energy c max ~ 10^ eV at the outer nonrelativistic shock created by the expanding 
cocoon around AGNs . 

The expected CR spectrum formed due their escape from the expanding shock has a 
form N oc e~ 7 with 7 ps 2.6 at energies e > 10 18 eV, which fulfills quite well the requirements 
for extragalactic CR sources. Therefore AGNs have to be considered as a prime candidate for 
the sources of extragalactic CRs. In such a case the observed CRs consist of two components: 
Galactic CRs, produced in SNRs, and extragalactic CRs from AGNs. 

The expected composition of ultra-high-energy CRs is characterised by two peaks in the 
energy dependence of the mean CR atomic number < A(e) >. The first one at the energy 



e ph 1 eV corresponds to the very end of the Galactic CR component (IBerezhko fc Volk 



20071 ). whereas the second, at the energy e ~ 10 19 eV, is expected at the beginning of the 
GZK cutoff. The strong energy dependence of CR composition within the energy interval 
from 10 17 to 10 18 eV is expected as a signature of the transition from Galactic to extragalactic 
CRs, whereas the detection of heavy CR composition at e ~ 10 19 eV has to be considered as 
the signature of CR production by nonrelativistic shocks. 
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